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1 Introduction 



The precise measurement of the V u t, element of the CKM quark mixing matrix from 
semileptonic b — > u decays is one of the most important goals for the B-factories. The 
high statistics accumulated by BaBar and Belle has recently made a measurement of the 
hadronic invariant mass distribution in inclusive B — > X u iv decays possible for the first 
time [1]. The new generation of analyses is based on fully reconstructed events, which 
allows high discrimination between charmless events and charmed background, even for 
hadronic invariant mass M x ^ 1.7 GeV. After unfolding detector and selection effects, 
BaBar has been able to measure the invariant mass distribution and two of its moments 
with promising accuracy. The measurements are possible even without an upper cut on 
M x , although it is clear that the relative error is smaller if one cuts at M^ ut close to 
the kinematic boundary for charm production (for instance, Ref.[l] adopts 1.86 GeV). 
The hope is that Mjp* can be raised enough to suppress non-perturbative effects that 
cannot be accounted for by the local Operator Product Expansion (OPE), namely Fermi 
motion effects related to the B meson distribution function(s), without compromising the 
experimental accuracy. Eventually, the measurement of the whole M x spectrum with this 
new experimental technique could provide complementary information on the distribution 
function and possibly a very clean extraction of V u b- 

Once accurate measurements of the moments of B — > X u lv distributions are available, 
the first task is to verify their consistency with moments of B — > X c lv and B — > X s 7 in 
the OPE framework. Present analyses of semileptonic and radiative moments [2] show an 
impressive consistency of all the available data and the non-perturbative parameters they 
provide agree with independent theoretical and experimental inputs. They determine the 
b quark mass within about 50 MeV and measure the expectation values of the dominant 
power suppressed operators with good accuracy. 

In this paper we extend a previous analysis of B — > X c tv moments in the kinetic 
scheme [3] to the B — > X u iv case. The moments in radiative decays in the kinetic scheme 
have been studied in [4]. Perturbative and non-perturbative corrections to the moments 
of B — > X u iv are given by the m c — > limit of those to the B — > XJLv moments (see 
[3, 5, 6] and Refs. therein), but the upper cut on the hadronic invariant mass, M™*, 
requires a dedicated study. We include all non-perturbative corrections through 0(1 /ml) 
[7, 8] and perturbative contributions through O(a 2 f3 ) [5, 6]. We also investigate the 
range of M^ ut for which the local OPE can be considered valid and give estimates of 
the residual theoretical uncertainty. A peculiarity of the b — > u case is the presence of 
a logarithmic divergence in the Wilson coefficient of the \jm\ correction [8, 9], related 
to the mixing between four-quark and the Darwin operators. Fortunately, the hadronic 
moments are relatively insensitive to the ensuing uncertainty. Finally, unlike the b — > c 
case, all perturbative and non-perturbative corrections to the b — > u case can be expressed 
in terms of simple analytic formulas, if only a lower cut on the electron energy is used. 

We also consider moments of the q 2 distribution. Experimentally, a measurement of 
the q 2 moments may soon become possible with only a lower cut on Eg. Provided E™ 1 
is not higher than say 1.5 GeV, the local OPE prediction should be reliable. That would 
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be a very interesting measurement, as the higher q 2 moments could efficiently isolate the 
effect of the Weak Annihilation (WA) contributions which are concentrated at high q 2 
values. From a practical point of view, as we will explain later on, the q 2 moments can 
be calculated using the same building blocks as the invariant hadronic mass moments. 

The paper is organized as follows: in Section 2 we briefly recall the formalism, de- 
fine our notation, discuss the kinematics involved in the experimental cuts, and describe 
the way we perform the calculation. We also provide tables of reference values and ap- 
proximate expressions, and discuss q 2 moments. In Section 3 we discuss the theoretical 
uncertainty of our results and the effect of Fermi motion. Section 4 summarizes our main 
results. In Appendix A and B we report analytic formulas for the non-perturbative and 
perturbative corrections to the moments in the case of a cut on the lepton energy only. 



2 The calculation 

We consider the normalized integer moments of the squared invariant mass, 

_ JdM 2 x M 2 x n dT/dM 2 x 

{Mx) ~ jdM 2 x dr/dM 2 x [i) 

and introduce the notation 

Ui = (M 2 x ), C/ 2 , 3 = ((M 2 X - (Mi)) 2 ' 3 ) (2) 

for the first three central moments. The physical hadronic invariant mass is related to 
parton level quantities by 

M 2 X = A 2 + 2m b AE + m 2 b s (3) 

where A = Mb — mi, (A is defined here to include all power suppressed terms), E Q = 
1 — v ■ q/rrib, and s = (v — q/rrib) 2 , g M is the four-momentum of the leptonic pair, and 
the four-velocity of the heavy meson. The moments of the parton level quantities E and 
s and of their product are obtained in the local OPE and are expressed in terms of the 
heavy quark parameters; in particular, they do not depend on M B or, equivalently, on A. 
The building blocks in the calculation of the moments in Eq. (2) are therefore 

1 r d 3 T 
M M = -j i E aiSa iE l siK 1E -^- IEt (4) 

where r = G 2 F m\ /(1927T 3 ) is the total tree-level width, (3 = 11 — 2/3 n/ with rif = 3, 
and M(jj) contains the tree-level contributions as well as non-perturbative corrections 
through 0(1 /ml). We compute the non-perturbative corrections using [7, 8], while the 
perturbative corrections A^K and are obtained (in the on-shell scheme) using the 
FORTRAN code accompanying Ref. [5] with a suitably small m c value. The one-loop 
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Figure 1: Effect of lepton energy and M x cuts on the E — s phase space. See the text for 
explanations. The cuts employed in the figure are Mjp* = 2.4 GeV and £ = 0.5. 

corrections A^-\ computed in this way agree with those computed from the results of 
Ref. [6], where the u quark is massless from the beginning. The numerical results for the 

(2) 

BLM corrections A\^^ are very sensitive to the value of the charm quark mass employed 

in the code [5], as at small m c the A$j\ are proportional to m 2 c In 2 m c . The numerical error 
associated with the choice of m c = 50 MeV in their computation is certainly acceptable 
for our purposes (it is below 1% in the BLM correction to the total rate, ^4( )> where it 
can be estimated using the exact result [10]). In the case where only a cut on the charged 
lepton energy is imposed, it is also possible to express M^j) and A^ in compact analytic 
form; the expressions relevant for the first three integer moments are given in Appendix 
A and B. In general, however, we rely on a numerical integration for the perturbative 
corrections. 

As mentioned in the Introduction, we consider here truncated moments subject to a 
lower cut on the energy of the charged lepton, E^ ut , and to an upper cut on the hadronic 
invariant mass, M™*. In the following we employ 

rpcnt 

* = 2^-. (5) 

m b 

It is useful to explain the kinematics with cuts in some detail. The region of integration 
in the E -sq plane is depicted in Fig. 1: the green (light) solid lines delimit the region 
of integration without any cut, that is the region between the curves Sq = 2E — 1 and 
s Q = Eq. The introduction of a cut in the lepton energy Ef nt divides this region into 
three parts that should be treated differently (see e.g. [11]): in the figure these regions 
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are separated by the black (dark) solid line that corresponds to s = (1 — £)(2E — 1 + £). 
In the first region, between s — 2-Eo — 1 an d s = (1 — £)(2E — 1 + £) (below the black 
line), one should use the differential rate calculated with the electron energy cut imposed. 
There are two regions above the black line, between s = (1 — 0(2-^0 — 1 + and s = Eq. 
In the lower of these two regions, the lepton energy is always above Ef^, and one should 
use the differential rate calculated without the cut. Finally, the upper region above the 
black line is excluded as the lepton energy is always below the cut. Whatever its value, a 
cut on the lepton energy affects both perturbative and non-perturbative contributions to 
the moments. 

A cut on the hadronic invariant mass M™* limits the region of integration to the area 
below the red dash-dotted line, that corresponds to M\ = (M™*) 2 . Increasing the value 
of M^ ut , the allowed region of integration expands. For a wide range of values of M^ 1 * 
and m&, as long as the red line does not come close to the E axis to the left of E — 1/2, 
the introduction of this cut affects only the perturbative corrections to the moments, as it 
excludes only events characterized by high s (hard gluon radiation). We limit ourselves 
to this case and consider only values of M™* above the lower limit 



corresponding to the situation in which the cut in M\ intersects the E axis at the value 
E = 1/2 (black dot). Using m b = 4.6 GeV, for instance, we have M£ ut > 1.89 GeV. 
In fact, the effect of the B distribution function becomes important within distances of 
O(Aqcd) from the E axis. As we will see in the next Section, a clean prediction of the 
moments requires Mjp* significantly higher than in the above equation. Table 1 reports 
our results for the components of M.{i,j) in the on-shell mass scheme at particular £ and 
M£ ut values. 

Although we start from on-shell expressions [5, 6], we actually employ the Wilsonian 
scheme with a hard factorization scale /i, whose optimal value is close to 1 GeV [12]. 
The Wilson coefficients in this scheme can be determined from the requirement that the 
observables be \x- independent. The initial condition is that at //— >0 one should recover the 
results in the on-shell scheme. In practice, at low perturbative orders this often reduces to 
re-expressing the pole-scheme results in terms of the running /i-dependent parameters. In 
particular, the /z-dependent parameters are the b quark mass mb(fj), kinetic expectation 
value /i 2 (/i), and the Darwin expectation value p| ) ( / u) 1 . In the following, our default choice 
for the non-perturbative parameters evaluated at \i = 1 GeV is 



The above values for the OPE parameters have been chosen having in mind the central 
values of the fits to B — > X c lv and B — > X s ^ moments [2] and some additional constraint. 

1 Unlikc Rcf. [3], we employ here a running Darwin expectation value. The relation between p^ D = p|,(0) 
and p D {li) can be found e.g. in [13]. 




(6) 



m b = 4.6GeV /4 = 0.40 GeV 2 fi 2 G = 0.35 GeV 2 
p^ = 0.1GeV 3 p| 5 = -0.1GeV 3 a s (m b ) = 0.22 . 



(7) 
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Table 1: Various contributions to the building blocks in the on-shell scheme with Ei > 1 GeV 
and M x < 2.5 GeV. The non-perturbative corrections in are calculated with the default 

values of the parameters given in Eqs. (7) and = 28. 



The coefficient function of the Darwin operator that contributes to the b^qlv total 



width at order \jm\ has a logarithmic divergence [8, 9] as m q ^0: 
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The singularity originates from the domain of low- momentum final-state quark (i.e., large 
and is removed by a one-loop penguin diagram that mixes the four-quark operator 
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^wa = 6&27o^l - ^l7°' u l — 2&27^l ^lH u< l (here in its Fierzed form) into the Darwin operator. 
Let us illustrate how this happens in the total semileptonic width: the lowest order 
contribution of the Weak Annihilation (WA) operator 0^ A is [16, 17] 



#r WA = r c WA {b\o^ a \b) 



(9) 



where C WA = 327i 2 /mf. In the factorization approximation the matrix element B WA = 
(B\0^ A \B) vanishes for O WA corresponding to zero lepton masses. Including 0(a s ) effects, 
the above equation becomes 



5F X 
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^WA 



1 + 0(a s )) B WA (fi 4q ) + ^ a (fi 4q )(B\bt a bJ2qiot a q\B) + ... 



CwA B WA (/i4g ) 
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(10) 



where \i 4q is the renormalization scale of the WA operator and a(// 4? ) is the contribution 
of a penguin mixing diagram renormalized in the MS scheme. We have used the fact that 
the Darwin operator is proportional to glbt a b Y, q Q 7o t a q by QCD equations of motion, 
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2.047 


u 2 


1.724 
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2.062 


2.228 


2.357 


2.377 


U 3 


1.188 


1.730 


2.338 


3.225 


4.198 


4.416 



Table 2: Hadronic moments for Ef nt = 1 GeV, B WA (0.S GeV) = 0, and the default values 
given in Eq. (7) at different M™*. The results are in GeV to the appropriate power. 

and we have neglected those contributions proportional to the matrix elements of 0^ A and 
of other operators that come with 0(a s ) Wilson coefficients, as they are irrelevant to the 
present discussion. The constant accompanying the logarithm in Eq. (10) depends on the 
renormalization scheme; it vanishes in the MS scheme that we have employed above 2 . We 
have therefore seen that the inclusion of the WA operator effectively replaces lnm^/m^ 
in Cd by lnm 2 //^ plus a constant, 

C D = - T -\Un4-^\. (11) 

Varying the renormalization scale /i^ adds a piece proportional to (? D to the WA expec- 
tation value. This contribution is independent of the flavor of the spectator, though, and 
therefore does not affect the differences between B + and B°. 

In the following we assume factorization to hold at the scale [i^ q = 0.8 GeV, i.e. 
_B WA (0.8GeV) ~ 0. A change in /i 4(J sets the natural size of the non-factorizable con- 
tribution in B WA . To get a crude estimate of how the non-valence (flavor-singlet) non- 
factorizable component of the expectation value of the WA operator affects the OPE 
predictions, we may vary /i^ in the interval 0.4GeV^/X4 9 ^1.7 GeV. It is clear, however, 
that ultimately the size of the WA expectation values in both B° and B + must be deter- 
mined experimentally. Notice also that, unlike the total width, the parton level moments 
M(jj) are not affected by the WA-Darwin mixing for i,j^0 (see Appendix A). 

In the calculation of the moments we follow [3] closely. In particular, we consider 
A = Oil) in the Aqcd/ 1 ^?) expansion and expand in Agc^/mfc and a s , neglecting all terms 
of 0(a s AQ CD /ml). This choice makes the hadronic moments sensitive to the choice of [x± q 

2 

in Eq. (11). One can parameterize this dependence using = 8 In ^ 2 L , and vary it in the 

range 16 < < 40, corresponding to the range in fi4 q just discussed, or equivalently use 
|i? WA (0.8 GeV) I < 0.004 GeV 3 . The contribution of the WA operator to the moments can 
be easily recovered in the following by the replacement p 3 D X^ — > p|, X^ — 32tt 2 B WA (ii4 q ). 

In Table 2 we provide some reference numbers for C/i )2 ,3, obtained using Ef nt = 1 GeV, 
X M = 28, and the default values given in Eq. (7) at different values of M™*. 

2 This applies if 0^, A is expressed in its Fierzed form in the continuation to D ^ 4 dimensions, which 
is part of the choice of scheme. Had we employed 0^ A = —ibljul u^b^ directly in D dimensions, 
using an anticommuting 75 (NDR scheme) , the logarithm In ml / n\ q would be accompanied by a constant 
+2/3. 
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Table 3: Coefficients of the linearized formulas in Eq. (12) for U\ at different cuts on the 
lepton energy without M x cut. 
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Table 4: Same as in Table 3 but for U 2 . 



In general, the BLM corrections are almost as relevant and have the same sign as the 
one-loop perturbative contributions at fixed a s = 0.22, i.e. they significantly decrease the 
effective scale of the QCD coupling. It is also convenient to have approximate linearized 
formulas for a generic moment M. of the form 

M(m b , /4, /4, pi, pI s - a s ) = V + B (m 6 -4.6 GeV) + A (a s - 0.22) 

+ P(^_o.4GeV 2 ) + D(p 3 D -0.1 GeV 3 ) (12) 
+ G (/4-0.35 GeV 2 ) + L (p 3 LS + 0.1 GeV 3 ) 
+ V5 WA (0.8 GeV); 

The values of V are obtained with the default values of the heavy quark parameters, 
B WA (0.8 GeV) = 0, and are quoted in GeV to the corresponding power. In Tables 3-5 we 
report values of the various coefficients with different E^ ut in the case without an M x cut. 
For values of M™* satisfying the bound in Eq. (6) only the perturbative contributions 
differ from the case without an Mx cut. Table 6 therefore shows only the coefficients V, 
B, and S for M^ ut = 2.5 GeV and Ej ut = 0.9 GeV. The results of our calculation are 
implemented in a FORTRAN code, available from the authors, that computes hadronic 
moments in b — > u for arbitrary Ef^ and for M™ 1 satisfying Eq. (6). 

We can compare the U\ values given in Tables 3 and 6 to the preliminary BaBar results 
of Ref. [1]. With a high M x cut of 5 GeV BaBar find lq(5 GeV) = 2.78 ± 0.82 GeV 2 , 
that is in good agreement with our reference value of 2.18 GeV 2 . BaBar also reports a 
result at low = 1.86 GeV which is very close to the lower bound of Eq. (6). In that 
case their result f/i (1.86 GeV) = 1.98 ± 0.20 GeV 2 is compatible with the reference value 
1.49 GeV 2 that we obtain at = 1.9 GeV, although Fermi motion effects, that shift 
Ui to higher values, have not been included in the calculation (see next Section). 
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Table 5: Same as in Table 3 but for U 3 . 
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Table 6: Coefficients of the linearized formulas in Eq. (12) for t/12,3 for M™* = 2.5 GeV and 
E^ ut = 0.9 GeV. Only the coefficients that differ from Tables 3-5 are reported. 



Finally, using our building blocks it is straightforward to study also q 2 moments. 
Indeed, replacing A with -mj in the rhs of Eq. (3) one obtains q 2 and can then calculate 
the q 2 moments 

^ _ J dq 2 q 2n dT/d q 2 
J dq 2 dT/dq 2 

in a way similar to the invariant mass moments. Table 7 gives the reference values and 
the coefficients of the linearized formula of Eq. (12) for the first three moments in the 
case of Ej ut = 1.2 GeV and no cut on M x . 



3 Theoretical uncertainty 

Let us now consider the various sources of theoretical uncertainty that affect our predic- 
tions. First we consider the uncertainty that affects the moments when no upper cut on 
M x is imposed. If the E e cut is not too severe (less than, say, 1.4 GeV), there are four 
main theoretical systematics: 

i) uncalculated 0(a 2 ) and O^Ag^/m^ 3 ) perturbative contributions to the Wilson 
coefficients; 

ii) missing 0(l/mf) non-perturbative effects; 

iii) the error from the scale in X M ; 

iv) Weak annihilation (WA) contributions. 

The first two items are common with the b — > c moments and can be analyzed in a similar 
way. The last two items, as clarified in the previous Section, are two facets of the same 
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Table 7: Coefficients of the linearized formulas in Eq. (12) for (g 2 ' 4 - 6 ) with E^ ut = 1.2 GeV. 

effect, and should not be counted twice. When we include the WA effects as a priori 
unknown, we use the variation of the scale in iii) to estimate the size of its flavor 
singlet contributions (flavor non-singlet WA effects can be studied from the difference in 
the moments for charged and neutral B). On the other hand, as discussed below, moment 
measurements allow to place constraints or to detect WA. In this approach fixing in 

extracts the expectation value of the WA operator normalized (in dimensional regu- 
larization) at this point. The remaining uncertainty comes from higher-order corrections 
to the Wilson coefficients, item i). 

In general, we note that b — > u moments are affected by larger theory errors than the 
b — > c moments. Using the simple recipe given in [3], we estimate the uncertainty related 
to i) and ii) above by varying p 2 , Pq, and a s by ±20%, p 3 D and p\ s by ±30%, and mj by 
±20 MeV in the theoretical predictions in an uncorrelated way. The typical results 

SU^Ui « 8% 5U 2 /U 2 « 20% 5U 3 /U 3 « 20%, (14) 

roughly reflect the theory errors due to i) and ii) above, independently of the accuracy with 
which we know the OPE parameters from [2]. They are driven by the strong sensitivity 
of Ui and U% to p 2 and p 3 D . 

The uncertainty of Eq. (14) can be estimated in alternative ways. For instance, we 
can evaluate C/1,2,3 using a different rearrangement of non-perturbative corrections, namely 
considering A as an O(Aqcd) quantity in the expansion in inverse powers of the b mass. 
In this case the moments are insensitive to X M and to the related error. The results are 
always within the ranges in Eq. (14). The main step necessary to improve on the above 
uncertainties is the calculation of perturbative contributions to the Wilson coefficients, of 
0(a 2 ) 3 and 0(a s /m 2 b ' 3 ). 

For what concerns the value of B WA , as mentioned in the previous section we vary it 
in the range |_B WA (0.8 GeV)| < 0.004 GeV 3 . This is a rather conservative estimate for the 
a priori unknown flavor singlet WA contribution, that induces a typical uncertainty of 
about 2%, 3%, 2% for Lq, U 2 , U 3 , respectively, although the error can be larger for high 
Ei cuts, as it is evident from Tables 3-5. 

We recall that WA contributions are concentrated at maximal q 2 , namely at the origin 
in Fig. 1. In the (E ,so) plane the fixed q 2 contours are straight lines identified by 
So = 2i?o — 1 + q 2 /m 2 . For instance, the r.h.s. boundary of the relevant phase space (the 
straight green line) corresponds to q 2 = 0. WA contributions are therefore characterized 

3 They are already available for the q 2 spectrum and moments [14], as well as for the total rate (first 
paper of [10]). 
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by small hadronic invariant mass and can be relevant in the total rate, but are suppressed 
in the moments of M\ . 

Turning to the theoretical uncertainty introduced by a cut on M x , we have already 
mentioned in the previous section that low M^ ut make the moments sensitive to the B 
distribution function. A rough but simple way to understand the range of M^ ut for which 
these non-perturbative effects become important is to rewrite Eq. (6) shifting A by the 
typical width of the distribution function, i.e. fi 2 /3 ~ ±0.37 GeV, to account for the 
Fermi motion. The result is that distribution function effects become important when 
M™ 1 is less than about 2.35 GeV. This result, however, is unlikely to apply to higher 
moments that are more sensitive to the tail of the distribution function. A detailed 
estimate would imply a dedicated implementation of the distribution function, which is 
beyond the scope of the present publication. 

A detailed study of similar effects in the photon energy moments of B — > X s + 7 
[4] showed that if the Wilsonian-type OPE with the hard cutoff at the scale around 
1 GeV is used, the inclusion of perturbative corrections affects only marginally the bias 
induced by Fermi motion. At the same time, the estimates are sensitive to the l/m,b 
corrections that decrease the variance of the light-cone distribution with respect to the 
heavy quark limit. Based on that experience, we have estimated the Fermi motion effects 
introduced by the Mx cut in the b^uiv moments by smearing the tree- level differential 
rate with an exponential distribution function (see eq.(13) of Ref. [4]) characterized by 
the low value /i 2 0.36 GeV 2 , to approximately account for the contributions to 

the second moment of the distribution function. Since the calculation of the hadronic 
moments is sensitive to the tail of the distribution, this estimate depends critically on 
the functional form adopted. In this respect, our choice of the exponential form leads to 
more conservative estimates than, say, with a Gaussian ansatz. Comparing the moments 
at M^ ut = Mb to those at various M™*, we find that the Fermi motion may alter Ui 
by ~ 6% at M^ ut = 2.5 GeV and by ~ 3% at M£ ut = 2.7 GeV. The higher moments 
are of course more sensitive: U2 may vary by 30% at Mjp* = 2.5 GeV and by ~ 20% at 
M£ ut = 2.7 GeV, while U 3 is dramatically affected - 0(1) effects - below M^ ut =2.7 GeV; 
the Fermi motion effects might become comparable to the other uncertainties only above 
Mf' = 3 GeV. Such high M^ ut values are certainly challenging for present experiments, 
but preliminary studies indicate the possibility of measuring the moments with interesting 
accuracy even for M™* well above the charm threshold [15]. 

Finally, let us discuss the uncertainty in the evaluation of the q 2 moments. We have 
already mentioned that these moments are very sensitive to the WA contributions, and 
therefore can be used to detect or place constraints on them. The flavor non-singlet WA 
contributions can be studied by comparing the q 2 moments (or other decay characteristics) 
measured in charged and neutral B decays separately [16]. For precision studies it is 
important to control the flavor singlet component of WA as well. They can be constrained 
by comparing the OPE predictions with data. In that respect, we need to consider only the 
uncertainties listed under i) and ii); regarding them in the same way as for the hadronic 
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moments, the results are 

S(q 2 )/(q 2 ) ^3% S(q 4 )/(q 4 ) w 6% S(q 6 )/(q 6 ) « 9% (15) 

and are dominated by perturbative effects. As mentioned in the previous section, the 
sensitivity of the various moments to the WA contributions can be understood from 
their dependence: varying £> WA (0.8 GeV) in the usual range and using the linearized 
formula of Eq. (12), we obtain a shift of approximately 4%, 11%, 21%, for the first three 
q 2 moments, respectively. Since we only consider q 2 moments without cuts on Mx, we do 
not have Fermi motion effects. 

4 Summary 

We have calculated the first three moments of the hadronic invariant mass distribution 
in charmless semileptonic decays in a Wilsonian scheme characterized by a hard cutoff 
/i ~ 1 GeV. Our calculation includes all known perturbative and non-perturbative effects, 
through 0(a 2 Po) and 0{l/m\) and is implemented in a FORTRAN code available from 
the authors. As required by the present experimental situation, we have considered cuts 
on the lepton energy and on the invariant hadronic mass and have obtained approximate 
formulas that summarize the dependence of the moments on the OPE parameters. 

The theoretical uncertainty of our OPE predictions ranges from 5% to 30%, but an 
upper cut on Mx introduces a dependence on the Fermi motion of the b quark in the 
meson. While we have performed a first estimate of these effects, the subject requires a 
more detailed investigation that we postpone to a future publication. Moreover, as the 
constraints on the shape of the distribution function are likely to improve in the future, 
our estimates of the Fermi motion effects should not be considered as an irreducible 
uncertainty. Conversely, the Mx spectrum and its truncated moments can themselves be 
used to constrain the distribution function, especially in the tail that is not accessible in 
radiative decays. 

We find that the bias introduced by the distribution function is not important for 
the first hadronic moment, if the M x cut is placed above 2.5-2.7 GeV. In that case its 
theoretical uncertainty is in the 10% range. The prediction for the second central moment, 
U2, is subject to a 20% uncertainty even without a cut on Mx- For cuts on Mx higher 
than ~ 2.7 GeV, the Fermi motion uncertainty on U2 may be as high as 20%. Finally, the 
third central moment, U3, is very sensitive to Fermi motion and can be predicted with a 
meaningful accuracy (~ 20%) only if the cut on M x is higher than ~ 3.2 GeV. 

We have also considered q 2 moments which could soon be measurable and are par- 
ticularly sensitive to WA. They can be predicted with good accuracy in the local OPE, 
provided the cut on the lepton energy is sufficiently low, E^ nt ^1.5 GeV. We have shown 
how it is possible to constrain both the flavor singlet and the flavor non-singlet WA con- 
tributions. In view of their potential interest, they should therefore be measured with 
Ep nt as low as possible. 
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Appendix A: non-perturbative corrections 

Here we give explicit expressions for the lowest order and OPE contributions to My, in 
the case of a lower cut on the charged lepton energy. They agree with Ref. [18] for £ = 0. 
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Appendix B: 0(a s ) corrections to the moments 



We report here analytic formulas for the 0(a s ) perturbative corrections to the building 
blocks My defined in Eq. (5) when a lower cut on the lepton energy is applied. The 
expressions are valid in the on-shell scheme for the b quark mass and they agree with 
Ref. [18] for £ = 0. We employ the short-hand = ln(l — £). 
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